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We study the spectral property of the antiferromagnetic Lipkin-Meshkov-Glick model. We show 
that at the supersymmetric point the spectrum of the model is gapped and its lower bound does not 
depend on the number of spins in the system. This implies that it is possible to generate specific 
entangled many-particle states in an ensemble of spins by adiabatic ground-state transitions in a 
scalable way. 


INTRODUCTION 


Fifty years ago, Lipkin, Meshkov and Click (LMG) proposed in a series of papers an exactly solvable model of 
interacting fermions in nuclear physics [ij. In recent years, due to the high controllability of trapped cold atomic and 
ionic systems, the model found applications in different areas of physics such as Bose-Einstein condensates Q, ion 
traps , in cavities Q . The model has also been used to investigate quantum phase transitions Q and their relations 
with bipartite entanglement ,1:0- Despite the fact that it has been proved that the model is integrable through 
the Bethe ansatz method 0 and some properties of the system can be extracted, in principle, this method does 
not yield an essential simplification for calculating the spectral properties of large systems. Among the theoretical 
methods applied to the LMG model were the WKB method [13 , the instanton technique [Il| , and a mapping onto a 
Schrodinger particle problem 0 . These methods are however only applicable in the asymptotic limit of large numbers 
of particles. 

For many applications, e.g. the creation of many-particle entanglement by adiabatic ground-state transitions, it is 
important to establish a lower bound of the gap between the ground state and the rest of the spectrum. For the finite 
number of spins it has been shown i, 0 that the antiferromagnetic supersymmetric LMG model ( ground state with 
zero energy and doubly degenerate excited states in the energy spectrum) is gapped for an arbitrary even number of 
particles. 

In this paper we present more detailed quantitative studies of the energy spectrum of this model. In particular, 
we find a lower bound for the spectral gap which does not depend on the number of particles. The latter result 
might sound paradoxical, because, it is commonly believed that the gap in the spectrum of spin-1/2 Hamiltonians 
decreases with increasing the number of particles. For local Hamiltonians in one spatial dimension with spin-1/2 it is 
an exact statement of the Lieb-Schultz-Mattis theorem 12|. In our model, however, the locality condition is violated 
and therefore this theorem is not applicable. 


COLLECTIVE SPIN COUPLING AND SUPERSYMMETRIC SPECTRUM 

We begin considering the LMG model. Let us consider an even number N of spin-1/2 particles interacting through 

N 

E (0 ("i) 

an , where cr/' denotes the i th component of the single-particle 

n—1 

spin. The interaction is assumed to be of second order in the total spin. The LMG Hamiltonian takes the following 
form 


H = ^ {xVl + XPI + AxiX2^) • 


( 1 ) 


For the interested reader about the LMG model we refer to the review article [13|. Hereafter, we consider only 
the maximum symmetric sector of the total Hilbert space, i.e. the total spin is equal N/2. The parameters Xi, 2 ,A 
are assumed to be positive constants. At the point where X 2 = 0, depending on the sign of the ^ the ground 
state is ferromagnetic or antiferromagnetic for ^ < 0 and ^ > 0 respectively. In the following we consider only the 
antiferromagnetic case. We note that by the rotation exp (i§ J^) the Hamiltonian m can be transformed into the 
form + XxiX2Jx and therefore, we may assume without loss of generality that xi > X2- As has been 

shown in [3| the spectrum of the antiferromagnetic LMG model becomes supersymmetric at the point A = 1. 
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FIG. 1: (color online) Typical supersymmetric spectra of the LMG model for J = 2 and J = 3 depending on the parameter 7 . 
The ground state zero energy not shown 


The analysis of the spectrum of the Hamiltonian operator (ED can be greatly simplified by introducing new param¬ 
eters, 

Xi = HoCOsh7, (2) 

X 2 = fio sinhy, 

where 

^o = xi-xl. tanh 7 =— (3) 

Xi 

and making use of the expression for a hyperbolic rotation about the x -axis with the parameter 7 

exp (-7 Jx) Jz exp (7 Jx) = Jz cosh 7 -|- iJy sinhy. (4) 

Thanks to this transformation, the Hamiltonian can be factorized as follows 

iT = fig (J 2 coshy -I- iJy sinhy) {J^ coshy — iJy sinhy) = 


= fig exp (-7 Jx) Jz exp (27 Jx) Jz exp (-yJx) 


(5) 


We set fig = 1 throughout. We see that for an arbitrary value of 7 all eigenvalues E of the Hamiltonian dS]) satisfy 
E > 0, and in addition to that for even number N of spin-1/2, i.e. integer angular momentum J = N/2, there is 
ground state with E = 0 : 


l^'o) 


1 

a/Pj (cosh 27 ) 


exp (7 Jx) \mz 


0 ) 


( 6 ) 


with Pj being Legendre polynomials (for more details, see Q and 0 ). 

Figured] shows the typical spectra for J = 2 and J = 3 depending on anisotropy parameter 7 . As one can see from 
Eq. ([Sj) the spectrum of the Hamiltonian is symmetric with respect to the 7 . 

A glance at the spectrum of the Hamiltonian as a function of 7 reveals (see FiglT|) that the spectrum consists 
of pairwise degenerate levels and a nondegenerate single state (the ground state) which signals the presence of a 
supersymmetry. In addition to that, we observe (Figl2|) that the gap does not close for large values N. The Hamiltonian 
([5]) has, for 7 7 ^: 0, no obvious symmetries other than the rotational symmetry around the a:-axis by ± 7 r and the spin-flip 
symmetry. It is clear that those symmetries cannot be responsible for this kind of spectrum. 

Now let us move on to an explanation of these observations. We begin by the Schwinger-boson representation of 
the angular momentum operators in terms of two bosonic operators a and b 


Jx = 

Jy = 

Jz = 2 ~ h%). 


( 7 ) 
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FIG. 2: (color online) Gap of the antiferromagnetic LMG model for different number of particles (J = 5,10,15, 25 and 30) as 
a function of anisotropy parameter 7 . The green line shows the analytical estimation (cosh 27 ) for the gap. 


For our purposes it is convenient to use another unitarily equivalent form of the Hamiltonian ([T|), namely 

H = cosh^ 1 + Jy sinh^ 7 + Jz cosh 7 sinh 7 . 

This Hamiltonian can be transformed into the following block diagonal form 


H = 


H 


even 


0 


0 

HodfZ 


( 8 ) 

(9) 


by grouping together its elements with even and odd bosomic excitations. In accordance with our partition Heven is 
a matrix with the dimension (n + 1) x (n + 1). It acts on states with even number of bosons such as (o’^) (b^) |0) 

where k is an even number. And Hom has n x n dimension and acts on states with odd number of bosons such as 
(a^)” (b^) ” |0) where n is an odd number. The Hamiltonian (| 8 ]) for iV = 2n particles has dimension {N + 1) x 

(TV + 1). The SUSY spectrum of the system must be twofold degenerate except for ground state i.e. the total number 
of states must be an odd number. Therefore we can see immediately that for an odd number of particles the SUSY 
should be broken. The decomposition @ allows us to introduce a ’’fermionic” number F: if F = 0 one has even 
excitations and F = 1 odd excitations. At the SUSY point (A = 1) one can introduce the following two off-diagonal 
(in the basis of odd and even bosonic excitations) operators (supercharges) 


Qi = o'xJx coshy -f cTyJy sinh 7 = 



(rt+l)x(r!,+ l) 

Jx cosh ^ + iJy sinh 7 


Jx cosh 'y — iJy sinh 7 



nxn 


( 10 ) 


Q 2 = cTyJx coshy — axJy sinh 7 = 


— Jy sinh 7 — fJa; cosh 7 

t+l)x(rt+l) 


— Jy sinh 7 -f i Ja; cosh 7 


0 


( 11 ) 


(n+l)x(n+l) 

One can check easily that the LMG Hamiltonian at the SUSY point can be represented as square of Qi ^2 namely 


H — Qi — Q2, 

and in addition to that Qi and Q 2 anticommute with each other 

{Qi, Q 2 } = 0. 

Combining relations (US and (TT^ in the form [l^ 

{Qi,Qj} = 25ijH, i,j = 1,2 
[Q^,H]=0, 


( 12 ) 

(13) 


( 14 ) 
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we arrive at a simple superalgebra, i.e., an algebra which includes both commutation and anticommutation relations. 

The above discussion can be made more concrete by looking at the case when iV = 4 (J = 2). In this case, the 
LMG Hamiltonian ([5|) after the renumbering states with even and odd excitations takes the following block diagonal 
form 


H = 


0 

0 0 

0 0 

The corresponding supercharges are 


= -27 


I 3 cosh 27 


3 

2 

=27 

0 

0 


0 

0 

0 

i (5 cosh 27 - 
3 


■ sinh 27 ) 


i (5 cosh 27 — sinh 27 ) 


Qi — 


0 

0 

0 

=-7 


0 

0 

0 

SpT' 


3g-7 


0 

0 

0 

0 

pT- 




3g-7 


and 


Q 2 = i 


0 

0 

0 

=-7 




^p-7 


-p-y 


0 


/3g7 /3„-7 


2 ^ 
-e'f 

0 

0 


A simple matrix calculation shows that these matrices satisfy the SUSY algebra (IT4l) . 

Because the Hamiltonian is Qf = Q\ its eigenvalues cannot be negative. Any states with zero energies are therefore 
ground states, and must be annihilated by both Q\ and Q 2 .- For the excited states, however, in gene ral, no further 
information can be gotten from the SUSY algebra ([T4|) . besides that they are two-fold degenerate [14 1. 

In the following we show that by using the representation ([5]) the degeneracy and an estimation for the gap can be 
obtained through direct calculations. 


SPECTRAL GAP OF THE LMG HAMILTONIAN 

It is easy to see from the Eq. ([5]) that the LMG Hamiltonian is similar to the non-hermitian Hamiltonian 
iJn = exp (-7 Ja:) H exp (7 Jx) = exp (-27 Jx) Jz exp (27 Jx) Jz = 

= JI cosh 27 -I- iJyJz sinh 27 , (15) 

From this form, one can see clearly that this Hamiltonian has a null state |mz =0) { Jz \mz = 0) = 0 for an integer 
angular momentum). The second term gives rise to a coupling between |toz) and \mz ± 1) states. The operator 
can be represented as a block matrix in a basis where Jz is diagonal, 

/iL_ 0 0 \ 

Hn = (a| 0 (a| , (16) 

Vo 0 H+J 

where iL+ and iL_ are permutation equivalent real J x J matrices, i.e. they have the same spectrum, (a| is the 
transposed vector which connects subspace with negative and positive values of the magnetic quantum numbers. For 
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our purposes, the generic form of (a| is not relevant. For example, the Hamiltonian for J = 2, takes the following 
form 


where 


and 




/ 4 cosh 27 sinh 27 0 0 0 \ 

— 2 sinh 27 cosh 27 0 0 0 

0 -^sinh27 0 -^sinh27 0 

0 0 0 cosh 27 — 2 sinh 27 

\ 0 0 0 sinh 27 4 cosh 27 / 


/ 4 cosh 27 sinh 27 \ 

— 2 sinh 27 cosh 27 J ’ ^ 


cosh 27 — 2 sinh 27 
sinh 27 4 cosh 27 


(a| 


-—sinh 27 , 



One can see that and are permutation equivalent, i.e. 


H_ = 




The matrix elements of with respect to the basis Jz are 


cosh ( 27 ) + ^ sinh 27 + “ ’^0 {J + m' + l)- + m') (J - m' + 1) 


(17) 


m,m' = —J,—J + 1,— 1. 


The spectrum of the can be obtained from the following algebraic equation 

Adet (A • Ijxj - H+) det (A • Ujxj - H_) = 0. 


(18) 


Since 77+ and i7_ have the same spectrum, the spectrum of is doubly degenerate except for the eigenstate 
jmz = 0). In other words, we have shown that the initial Hamiltonian ([5]) in the case of integer angular momentum 
has a supersymmetric spectrum Q . Hence, to find the excited spectrum of the LMG Hamiltonian it is enough to find 
the spectrum of the 77+ or 77_. The ground state energy of the 77+ or 77_ will corresponds to the gap of the original 
Hamiltonian. 

As one can see from the Eq dm, the matrix 77_ is a real tridiagonal matrix and its upper and lower sub-diagonals 
have opposite sign. 

Let us consider a real tridiagonal matrix with real spectrum 


ai —Pi 
7i a2 —P 2 


—Pn-i 

'Jn — l 


(19) 


with the matrix elements 


Pi,'Yi >0, * = 1, 2, ...n - 1. 


( 20 ) 


By the similarity transformation 


A' = TAT-^, 


( 21 ) 
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one can transform this matrix into a more symmetric form 


A' = 


V/5i71 —VP 2 I 2 

V/5272 


■ _ —\/ Pn-lln-l 


( 22 ) 


where T is a diagonal matrix with diagonal entries 

h = 1 , 

The real matrix A' can be written as 


U 

ti+i 


A' + A'^ A'-A'^ 

A' = -h i -, 

2 2 i ’ 


(23) 


(24) 


where the real part ^ ^ is diagonal. Let A„ {A') = A„ (A) > 0 and \tfn) are smallest eigenvalue and corresponding 

normalized eigenvector of the A'. An estimation for A„ can be obtained from the decomposition (IMl) . Indeed, 


A' _|_ J^'T 

Xn (A) = {ipn\A' \ipn) = {^n \ -^- Wn) > OLn = min(afc) 

Applying this estimation, we obtain the desired result a lower bound for the gap 

A A = En {H-) > cosh 27 

between the first two eigenvalues of the antiferromagnetic supersymmetric LMG model. 


(25) 

(26) 


DISCUSSION 

We have shown that the antiferromagnetic supersymmetric LMG model with even number of spins is gapped. 
Surprisingly the obtained lower bound of the gap does not depend on the total number of spins in the system. This 
means that in contrast to one dimensional local systems with spin- 1/2 where the gap decreases with increasing 
number of particles for spins-1/2 systems interacting through a collective coupling of the Lipkin-Meshkov-Glick type 
this is in general not the case. Figure [2] shows that our estimate is almost exact for small 7 . Therefore, one can think 
naively that the result (j26ll for small 7 is a perturbative estimate. The above consideration shows, however, that usual 
perturbation theory is not applicable at the supersymmetric point. We note that for small 7 the LMG model can be 
considered as an example of a spin tunnelling system 0 - In order to calculate the gap of the spectrum for small 7 , 
one has to find the tunnelling splitting of pairwise degenerate excited states. The supersymmetric models are special 
ones due to the fact that in any order of 7 the tunnelling splitting must disappear. It is worth noting, however, that 
even if it is would be possible to find some estimation for the small 7 (by some perturbative arguments), it will still be 
insufficient for applications such as creation of many-particle entangled states by adiabatic ground-state transitions, 
where an adiabatic process starts at 7 —> 00 ( ground state is a product state) and end up in an entangled state at 
7 = 0 . 

An interesting extension of this work will be to study the LMG model with broken supersymmetry (A ^ 1) and 
the dependence of the adiabatic condition on the size of the system. A detailed study of this problem will be given 
elsewhere. 

I am grateful to M. Fleischhauer for many fruitful and stimulating discussions. 
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